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ABSTRACT 

Following the suggestion of Black(1997) that some massive extrasolar planets may 
be associated with the tail of the distribution of stellar companions, we investigate a 
scenario in which 5 < N < 100 planetary mass objects are assumed to form rapidly 
through a fragmentation process occuring in a disc or protostellar envelope on a scale 
of 100 au. These are assumed to have formed rapidly enough through gravitational 
instability or fragmentation that their orbits can undergo dynamical relaxation on a 
timescale ~ 100 orbits. 

Under a wide range of initial conditions and assumptions the relaxation process 
ends with either (i) one potential 'hot Jupiter' plus up to two 'external' companions, 
i.e. planets orbiting near the outer edge of the initial distribution; (ii) one or two 
'external' planets or even none at all; (iii) one planet on an orbit with a semi-major 
axis 10 to a 100 times smaller than the outer boundary radius of the inital distribution 
together with an 'external' companion. Most of the other objects are ejected and 
could contribute to a population of free floating planets. Apart from the potential 
'hot Jupiters' , all the bound objects are on highly eccentric orbits. We found that, 
apart from the close orbiters, the probability of ending up with a planet orbiting at 
a given distance from the central star increases with the distance. This is because of 
the tendency of the relaxation process to lead to collisions with the central star. 

We discuss the application of these results to some of the more massive extrasolar 
planets. 

Key words: giant planet formation - extrasolar planets - dynamical relaxation - 
orbital elements 



1 INTRODUCTION 

The recent discovery of extrasolar giant planets orbiting 
around nearby solar-type stars (Marcy & Butler 1998, 2000; 
Mayor & Queloz 1995) has stimulated renewed interest in 
the theory of planet formation. The objects observed so far 
have masses, M p , that are characteristic of giant planets, i.e. 
0.4 Mj < M p < 11 Mj, Mj denoting a Jupiter mass. The or- 
bital semi-major axes are in the range 0.04 au < a < 2.5 au, 
and orbital eccentricities in the range 0.0 < e < 0.67 (Marcy 
& Butler 2000). 

It is a challenge to formation theories to explain the ob- 
served masses and orbital element distributions. There are 
two main theories of giant planet formation (see Papaloizou, 
Terquem & Nelson 1999 and references therein). One is the 
core instability scenario. In this, a solid core of several earth 
masses is built up in the protostellar disc, at which point it 
is able to begin to accrete gas and evolve to become a gi- 



ant planet. Once massive enough it is able to open a gap in 
the disc and undergo orbital migration through disc proto- 
planet interactions (e.g. Lin & Papaloizou 1993). It has been 
suggested that the high orbital eccentricities for extrasolar 
planets might be explained by disc-protoplanet interactions 
(Artymowicz 1992). 

Recent simulations of protoplanets in the observed mass 
range (Kley 1999, Bryden et al. 1999, Lubow, Seibert & 
Artymowicz 1999) interacting with a disc with parame- 
ters thought to be typical of protoplanetary discs, but con- 
strained to be in circular orbit, indicate gap formation and 
upper mass limit consistent with the observations. However, 
simulations by Nelson et al. (2000) that relaxed the assump- 
tion of fixed circular orbits found inward migration and that 
the disc-protoplanet interaction leads to strong eccentricity 
damping. Thus the observed eccentricities of apparently iso- 
lated extrasolar planets are so far unexplained by this sce- 
nario. 
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The other possible formation mechanism is through 
fragmentation or gravitationnal instability in a protostel- 
lar disc (e.g. Cameron 1978, Boss 2000). This may occur 
early in the life of a protostellar disc surrounding a class 
protostar on a dynamical timescale. Such discs have been 
observed (see, e.g., Pudritz et al. 1996) and the character- 
istic size is about 100 au. It is unlikely that such a process 
would operate at distances smaller than about 50 au from 
the central star as, in the optically thick parts of the disc, 
non axisymmetric density waves redistibute mass and angu- 
lar momentum before fragmentation can proceed ( eg. Pa- 
paloizou & Savonije 1991, Laughlin & Bodenheimer 1994). 
Fragmentation is more likely when cooling is efficient, as 
may occur in the optically thin parts of the disc, beyond 
about 50 au (Papaloizou et al 1999). However, the detailed 
conditions required for it to occur are unclear and may re- 
quire constraining influences from the external environment 
( Pickett et al 2000). Note that fragmentation may also oc- 
cur before a disc is completely formed, during the initial 
collapse of the protostellar envelope. Such opacity limited 
fragmentation has been estimated to produce objects with a 
lower mass limit of 7 Jupiter masses (Masunaga & Inutsuka 
1999), but there is no definitive argument to rule out some- 
what smaller masses (Bodenheimer et al. 2000) . It is possible 
that both a disc and fragments may form simultaneously out 
of the envelope, the relative importance of the two processes 
depending for instance on the angular momentum content of 
the envelope, on the strength of any magnetic field ( so far 
neglected in disc fragmentation calculations) and possibly 
on the initial dumpiness. Note that large scale observations 
of class envelopes so far do not rule out the presence of 
clumps with masses smaller than about 10 Jupiter masses 
(Motte & Andre 2001). 

It is the purpose of this paper to investigate the evo- 
lution under gravitational interactions of a distribution of 
N massive planets which we assume to have been formed 
through a fragmentation process rapidly enough that their 
orbits can undergo subsequent dynamical relaxation on a 
time scale of hundreds of orbits. In common with related 
work on orbital evolution ocuring after assummed formation 
in a disc (eg. Weidenschilling & Mazari 1996, Rasio & Ford 
1996, Lin & Ida 1997) we shall neglect the effects of any rem- 
nant disc gas so that apart from tidal interactions with the 
central star there are only gravitational interactions. Thus 
this work complements studies of the initial fragmentation 
process in a gaseous medium. 

It turns out that the resulting evolution leads to similar 
end states independently of whether the initial configuration 
is assumed to be in the form of a spherical shell or a disk-like 
structure. 

The motivations for this work are firstly the suggestion 
by Black (1997) and Stepinski & Black (2000) that massive 
extrasolar planets on highly eccentric orbits could actually 
be the low-mass tail of the low-mass companion distribu- 
tion to solar-like stars produced by fragmentation processes. 
Here we wish to investigate to what extent planets with or- 
bital elements similar to those observed can be produced 
and in particular whether 'hot Jupiters' orbiting close to 
the star can be formed. Secondly we consider the recently 
detected population of free-floating planets and its relation- 
ship to that of planets orbiting solar-type stars ( Lucas & 
Roche 2000, Zapatero Osorio et al. 2000). It is of interest to 



know to what extent free-floating planets could be produced 
as a result of ejection from the neighbourhood of a star. 

We have considered the orbital evolution of N bodies 
with masses in the giant planet range, which are assumed 
to be formed rapidly, using up the gas in a protostellar disc 
or envelope around a solar mass star, so that they can un- 
dergo subsequent dynamical relaxation on a timescale ~ 100 
orbits. We have performed calculations with 5 < N < 100. 

In all the runs we performed, most of the planets where 
ejected from the system and at most 3 planets remained 
bound to the central star. We found that close encounters 
with the central star occured for about 10% of the planets 
for all values of N considered. Such close encounters early 
in the evolution tended to result in collisions. These tended 
to be avoided at later times so that tidal interaction might 
then result in orbital circularization at fixed pericentre dis- 
tance leading to the formation of a very closely orbiting giant 
planet. 

Typically the runs ended up with either (i) one potential 
'hot Jupiter' plus up to two 'external' companions, i.e. plan- 
ets orbiting near the outer edge of the initial distribution; 
(ii) one or two 'external' planets or even none at all; (iii) 
one planet on an orbit with a semi-major axis 10 to a 100 
times smaller than the outer boundary radius of the inital 
distribution together with an 'external' companion. Apart 
from the potential 'hot Jupiters' , all these objects are on 
highly eccentric orbits. We found that, apart from the close 
or biters, the probability of ending up with a planet orbiting 
at a given distance from the central star increases with the 
distance. 

The objects that become unbound may contribute to a 
population of freely floating planets ( Lucas & Roche 2000, 
Zapatero Osorio et al. 2000) which could be several times 
larger than that of giant planets found close to the central 
star. 

Thus the dynamical relaxation process considered here 
may operate in some cases to produce giant planets with 
high orbital eccentricity at several astronomical units from 
their central star as well as closely orbiting planets. In all 
cases a population of loosely bound planets is also expected. 

The plan of this paper is as follows. In section [| we 
describe the model and basic equations used. In section ^ 
the initial conditions are formulated and in section k the 
physics of the relaxation process is discussed. In section 1 wc 
present our numerical results and in section ^| we summarize 
and discuss them. 



2 THE MODEL AND BASIC EQUATIONS 

We consider a system consisting of iV planets and a primary 
star moving under their gravitational attraction. 
As we are interested in possible close approaches to, or col- 
lisions with, the central star we adopt a spherical polar co- 
ordinate system (r, 9, <p) with origin at the stellar centre. 
The planets and central star are treated as point masses. 
However, to take into account possible losses of orbital en- 
ergy and angular momentum to the stellar material, a simple 
model for taking into account the tidal interaction between 
the star and a closely approaching planet is also included. 
The equations of motion are 
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dt 2 



GM,Vi \ - GMjVi 



1 



(1) 



Here M*, Mj, rj and ri,- denote the mass of the central star, 
the mass of planet i, the position vector of planet, i, and r< — 
rj respectively. The acceleration of the coordinate system 
based on the central star (indirect term) is 



\ - GM 3 Vj 

|p 3 -|3 



(2) 



and that due to tidal interaction with planet i is Ft*. 

In the situation envisaged here, tidal interactions occur 
when a planet has a close encounter with the star. When 
this occurs, the planet approaches from large distances on 
an almost parabolic orbit. The time between a subsequent 
encounters will then be long compared to that for the tidal 
interaction itself. Accordingly we approximate the process 
as a sequence of independent energy and angular momen- 
tum transfers that occur at each periastron passage. We 
utilize the results of Press & Teukolsky (1977) who calcu- 
lated these transfers in the small perturbation limit for a 
non rotating star modelled as a polytrope and a perturber 
on a parabolic orbit. We shall neglect the effects of tides 
acting on the planet itself. Accordingly our model is sim- 
plified. However, it does enable us to include tidal effects 
and demonstrate how they start to lead to orbital circular- 
ization and gravitational decoupling of an inner planet from 
the others as it moves onto a close orbit. However, it is only 
applicable while the planet orbit has high eccentricity. 

We adopt a form for Fxi that is able to approximately 
give the correct angular momentum and energy exchanges 
with the star on a close approach but which is negligible at 
larger distances from the central star. 



GMi R* Rpi T\ 



Ci\ji\\n 



i |2 

ri v. 



Vi ■ ViVi 



(3) 



Here ji is the specific angular momentum of planet i, 
R, is the stellar radius, R P i — jf /(2GM*) is the distance 
of closest approach of planet i assuming a parabolic orbit, 
CI = 20F/3, and Ti = 0.6/(1 + (R pl /R,) 3 ). 

The equations of motion are here integrated using the 
Bulirsch-Stoer method ( eg. Press et al 1993). 

Using we can derive the energy lost to the star dur- 
ing a close encounter of planet i, assumed on a parabolic 
orbit as 



AE = - 



(4) 



where, because of the rapid convergence of the integral, the 
limits are extended to ±oo. This gives 



AE = 



R 6 - 



(5) 



which gives values coinciding approximately with values 
given by Press & Teukolsky (1977). 

We comment that acording to (|H|) a star grazing en- 
counter on a parabolic orbit results in a final semi-major 
axis a ~ 1.7R»M»/Mi. For a central solar mass with R, = 
10 11 cm. and Mi ~ lMj, a ~ 10 au. Thus bound orbits 
with a ~ 10 au are significantly modified if they have a 



close encounter with the central star. Note too that the en- 
ergy exchange rates are small for the planetary mass objects 
considered here giving some justification for the linear ap- 
proximation used to calculate them. 



3 INITIAL CONDITIONS 

The simulations performed here were begun by placing N 
planets in some specified volume in a random location cho- 
sen to give a Monte-Carlo realization of a prescribed den- 
sity distribution. We considered both the case of a uni- 
form density spherical shell with R m in < r < Rmax and 
that of a thick annulus with R m in < r < Rmax occupying 
cos _1 (0.1) < 9 < cos~ 1 (— 0.1) with a density oc r -2 , where 
r is the spherical radius and 9 the co-latitude. For the cal- 
culations presented here, Rmin = O.lRmax and 0.5 R ma x 
for the annulus and the spherical shell, respectively. The 
planets were then given the local circular velocity in the 
azimuthal direction. Note that because of the spatial place- 
ments the initial orbits are not coplanar. In some cases the 
planets were taken to have equal mass M p while in others 
each planet was allocated a mass qM p where q was a random 
number between zero and one. In the latter case the mass 
distribution of the protoplanets does not correspond to the 
specified density distribution indicating some redistribution 
of mass among the embryos. 

In general the various initial conditions we have consid- 
ered lead to the same qualitative evolution. 

However, it must be emphasized that the systems 
are chaotic with the consequence that very small detailed 
changes to the initial conditions or the integration proce- 
dure will in general lead to very different results in detail. To 
deal with this issue one can adopt the notion of shadow or- 
bits (Quinlan & Tremaine 1992). According to this, one can 
expect that although inevitable small errors lead to a signifi- 
cant deviation of the numerical solution from the actual one, 
there is another solution of the real system obtained with 
slightly different initial conditions that remains close to the 
numerical one. Thus we should be able to identify qualita- 
tive trends in the evolution of an ensemble of systems of the 
type we consider and we restrict ourselves to that. 



3.1 Scale Invariance 

We comment that the equations we solve, incorporating tidal 
effects, have a radial scale invariance. Thus all radii may be 
mutiplied by some factor / while the timescales are multi- 
plied by / 3/2 . The size of the central star has to be scaled 
by the factor / also. Here we shall take the unit of length 
to be Rmax- Then the stellar radius is specified through 
R*/Rmax- The time unit is the period of a circular orbit 
at Rmax, Pq = 2-KR^lax / V GM* . The interactions amongst 
the planets lead to some escaping the system. Objects with 
positive energy which had reached a distance /3R m ax were 
considered to be escapers. We have considered (3 = 20 and 
/3 = 100 which both lead to the same qualitative picture. 
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4 EVOLUTION RELAXATION AND STAR 
GRAZERS 

Even the systems which have a small number of bodies are 
found to interact strongly and to undergo relaxation like 
a stellar system (Binney & Tremanine 1987). For such a 
system the realaxation time is 

°- 34 " 3 (6) 



tR 



3V3G 2 M P pln(A) 



Here the root mean square velocity is v, p is the mass density 
of interacting bodies assumed for simplicity to have equal 
mass M p , and A = M»/M p . 

Using v 2 = GM,/R, NM V = A-kR^p/Z and adopting 
the orbital period P — 2-kR/v this becomes 



t R = 



0.0UM 2 P 
M 2 N\n(M t /M p ) 



(7) 



Thus for M p /M, = 5.0 x 10" 3 and N = 5 the relaxation 
time is about 100 orbits. For systems with R in the 100 au 
to 1000 au range this time is around 10 5-6 y which is within 
the estimated lifetimes of protostellar discs. 

The evolution we obtain is similar to that undergone 
by spherical star clusters (Binney & Tremaine 1987). The 
relaxation, due to binary encounters, causes some objects 
to attain escape velocity while others become more bound. 
Eventually all either escape or end up in extended orbits, 
except for one body which takes up all the binding energy. 
This is a generic result provided that close encounters with 
the star can be avoided. However, for the parameter range 
of interest such encounters are likely. The situation resem- 
bles that of accretion of stars from a spherical star cluster 
by a black hole (eg. Frank & Rees 1976) . At a location with 
radius R, the fractional volume of phase space containing 
orbits that would collide with the star if unperturbed is 
~ R*/R ( this being the ratio of the square of the angu- 
lar momentum below which an impact is expected to the 
square of the mean angular momentum). If an object can- 
not diffuse out of this volume before impact, the mean time 
before diffusion into the effective volume produces an im- 
pact, for a particular object, is ~ tR. However, the time to 
diffuse out of the effective volume is ~ (R*tn.)/R. If this 
time is less than the crossing time t c ~ R/v, the expected 
mean time to impact is t enc = (Rt c )/R*. This time is just 
the crossing time divided by the probability of being in the 
effective volume of phase space which is regularly sampled 
because of the effective diffusion (see Frank & Rees 1976). 

In our calculations the innermost most tightly bound 
object undergoes relaxation or phase space diffusion which 
can lead to close encounters. For R = 25 au, t c ~ 20 y, and 
P* ~ 1.5 x 10 11 cm, t enc ~ 5 x 10 4 y. Thus we can ex- 
pect close encounters to occur within the general relaxation 
process. In some cases these can lead to a strong tidal inter- 
action which can circularize the orbit and potentially lead 
to the formation of a 'hot Jupiter'. 



5 NUMERICAL RESULTS 

We here describe a sample of our results which illustrate 
the characteristic behaviour exhibited by the systems we 
consider. The calculations presented are listed in table [ij 



Run 


N 


/ Rmax 


In 


M p /M, 




1 


5 


0.0 


12 


5 x 10~ 3 




2 


10 


9.396 x 10" s 


11 


5 x 10- 3 




3 


10 


1.337 x 10~ 4 


12 


5 x 10~ 3 




1 


8 


1.337 x 10~ 4 


12 


5 x 10~ 3 


(a) 


5 


10 


9.396 x 10- s 


11 


5 x 10~ 3 




6 


100 


9.396 x 10" s 


11 


1 x 10~ 2 


(R) 



Table 1. This table lists the parameters for each simulation: the 
number of planets TV, the ratio of stellar radius to Rmax and the 
initial set up In, where n = 1, 2 denotes spherical shell and thick 
annulus, respectively. The last column contains Mp/M*, with (R) 
denoting that the masses were allocated uniformly at random in 
the interval (0; M p ). 



We note that if Rmax =100 au, R,/R max = 9.396 x 10" 5 or 
1.337 x 10 -4 corresponds to P, — 2 or 3 Rq, respectively, 
which is the radius of a protostar with a mass around 1 Mq 
in the early stages. 

Run 1 corresponds to a system of N = 5 planets each 
having a mass 5 x 10 _3 M*. Figure |l| shows plots of a/R max 
in logarithmic scale, where a denotes the semi-major axis, 
for each planet in the system versus time (measured in units 
of Po). Each line corresponds to a different planet. During 
the run 3 planets escape, a line terminating just prior to an 
escape. For this case P* = so that there was no tidal inter- 
action with the central star. The initial relaxation time for 
this and other similar cases is on the order of 100Po , in agree- 
ment with equation (^). The evolution of a(l — e)/Rmax, 
a(l — e) being the pericentre distance and e being the eccen- 
tricity, is also shown in Figure [j]. We see that approaches to 
within ~ W~ 2 Rmax of the central star occur for the inner- 
most object on a timescale of the same order of magnitude 
as the relaxation time. This is a common feature of the sim- 
ulations presented here. Note that for Rmax = 100 au, the 
closest approach is to within ~ 3 x 10 11 cm, which is com- 
parable to a solar radius. After about 6, 000Po, the main 
relaxation is over with nearly all the binding energy being 
contained within one object with a ~ O.lRmax and e ~ 0.9. 
However, After about 5 x 10 4 Po, the two innermost planets 
have a close approach which results in their position being 
exchanged and one of them being ejected. The end result of 
this run is then two planets on well separated orbits. The 
innermost planet has a ~ O.lRmax and e ~ 0.5 

A plot of the evolution of the semi-major axes and peri- 
center distances versus time (measured in units of Po) for 
run 2 in table |l| is given in Figure ^. For this run the initial 
number of planets is N = 10 and R*/R m ax = 9.396 x 10~ J so 
that tidal effects operate. After about 10 4 Po, there are only 
two planets still bound to the star. The others have been 
ejected. The innermost planet which is left has e = 0.66 
and a ~ 0.02R ma x- If Rmax = 100 au, this corresponds to 
a ~ 2 au. 

Figure ^| shows the evolution of the semi-major axes 
and pericenter distances versus time (measured in units 
of Po) for run 3 in table [l]. For this run N = 10 and 
R*/Rmax = 1.337 x 10 -4 . The plot terminates at about 
3.4 x 10 4 Po, just before the innermost planet has a close 
encounter with the star. The evolution toward the end of 
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the run is zoomed on in Figure ]A We see that the semi- 
major axis decreases significantly, down to about 10 -2 Rmax, 
whereas the pericenter distance varies much less. This indi- 
cates a tidal interaction with the star. The run if contin- 
ued much further becomes innaccurate because of our crude 
treatment of tides, as mentionned in section ^. We expect 
at that point the orbit to circularize at fixed pericentre dis- 
tance. Evolution of this type which is also quite common in 
our simulations could potentially lead to a 'hot Jovian mass 
planet'. 

Figure ^ shows the evolution of the semi-major axes and 
pericenter distances versus time (measured in units of Po) 
for run 4 in table |l| In this case N = 8,R* / Rmax = 1.337 x 
10" 4 and the masses were selected uniformly at random in 
the interval (0; 5 x 10 -3 M*). Here again, after about 3 x 
10 4 P there are only two planets still bound to the star. 
The innermost planet has e = 0.25 and a ~ 0.04R m ax, which 
would correspond to 4 au in a 100 au annulus. 

We now consider runs with larger N. Figure ^| shows 
the evolution of semi-major axes and pericenter distances 
versus time (measured in units of Po) for run 5 in table |l| 
In this case N = 40 and R»/R ma x = 9.396 x 10" 5 . In that 
case also there are only two planets bound to the star after 
about 8, OOOPo- The others have either been ejected or have 
collided with the central star. Our treatment of tides is too 
crude to be certain about the orbital evolution of planets on 
very eccentric orbits which have grazing approaches to the 
star. In our simulations, such approaches result in the planet 
being lost. As far as the subsequent evolution of the system 
is concerned however, it does not matter whether the planet 
indeed hits the star of gets circularized on a close orbit. 
In this run, evolution occurs on a shorter timescale due to 
the larger number of planets. Here the innermost planet has 
e ~ 0.9 and a ~ 0.03R ma x, which is 3 au if Rmax = 100 au. 

Figure |7| shows the evolution of the semi-major axes 
versus time (measured in units of Po) for the only 2 planets 
in the system which have not either been ejected or col- 
lided with the central star after about 180 Po, for run 6 
in table [[]. Also shown are the pericenter distance and 
the apocenter distance for these planets. Here N = 100, 
R,/R max = 9.396 x 10 -5 and the masses were selected uni- 
formly at random in the interval (0; 10 _2 M«). We see from 
Figure |^ that at a time of about 1.5 x 10 4 Po the apocenter 
and the pericenter of the innermost and outermost planets, 
respectively, are very close to each other. The innermost 
planet then suffers a gravitational scattering by the more 
massive outermost planet which results in an increase of 
its eccentricity and eventually a collision with the central 
star. We are then left with one planet with e = 0.85 and 

O ~ 0.7 Rmax- 



6 SUMMARY AND DISCUSSION 

We have considered the orbital evolution of N bodies with 
masses in the giant planet range, which are assumed to have 
formed rapidly enough that they can undergo subsequent 
dynamical relaxation on a timescale ~ 100 orbits. 

We have considered 5 < TV < 100. We assume that rapid 
formation may occur as a result of fragmentation triggered 
by gravitational instabilities or clumping, either in a spher- 
ical envelope during the initial protostellar collapse phase 



or in a disc like configuration (e.g. Pudritz et al 1996) as 
that forms. But note that as magnetic fields may play a role 
the disc may not be entirely centrifugally supported. We 
have considered initial distributions of planets with masses 
in the range 5 to 10 Jupiter masses that are in the form of 
a sphericalshell or are disc like. However, final outcomes are 
independent of this. In our calculations, along with previous 
work related to the core accretion model (Weidenschilling & 
Mazari 1996, Rasio & Ford 1996, Lin & Ida 1997) we ne- 
glect the effects of gas, thus the efficiency of fragmentation 
is assumed to be maximal. 

Although our results can be applied to different radial 
scales, we focus the discussion below on distributions with 
an outer radius of 100 au this is the dimension of observed 
protostellar discs in class objects (eg. Pudritz et al 1996). 

This work has been partly motivated by the sugges- 
tion by Black (1997) and Stepinski & Black (2000) that 
some massive extrasolar planets could actually belong to the 
low-mass tail of the distribution of low-mass companions to 
solar-like stars. That suggestion is derived from the obser- 
vation that extrasolar planets far enough from the star that 
tidal circularization does not operate tend to have highly 
eccentric orbits and that these could be produced by a re- 
laxation process of the type we consider. 

Altogether we have run about 25 cases with 30 < TV < 
100 and many more with N < 10. We have described six 
representative runs in detail in the paper. Independently of 
how many planets we began with, the initial set up and the 
technical details of the computation, we obtained the same 
sets of characteristic behaviour and end states. In every case 
most of the planets where ejected from the system and at 
most 3 planets remained bound to the central star after a 
time typically on the order of a few 10 4 outer periods Po. 
The dynamical relaxation phase is shorter when the initial 
number of planets is larger, but when the number of planets 
is reduced to less than 10 the system evolves in the same 
way as the systems which had a smaller intital number of 
planets. 

We found that close encounters with the central star 
often occured (for about 10% of the planets) for all values 
of N considered. At an early stage these tended to result 
in direct collisions. When a direct collision is avoided, tidal 
interaction between the star and planet on a very ecentric 
orbit may result in orbital circularization at fixed pericentre 
distance which might ultimately lead to the formation of a 
very closely orbiting giant planet. 

As far as we could monitor the runs presented here, 
we did not find any physical collisions between the planets 
themselves. This is consistent with Lin & Ida (1997) who 
found that such collisions were very rare when, as here, mu- 
tual tidal interactions between the planets were assumed 
to be ineffective. If physical collisions were to occur , they 
would not be expected to affect the typical outcome of our 
runs. 

Typically the runs ended up with either (i) one poten- 
tial 'hot Jupiter' plus up to two 'external' companions, i.e. 
planets orbiting near the edge of the initial distribution; (ii) 
one or two 'external' planets or even none at all since a grav- 
itational scattering between two planets may result in one 
being ejected and the other one colliding with the star; (iii) 
one planet on an orbit with a semi-major axis 10 to a 100 
times smaller than the inital distribution, e.g. a ~ 1-10 au 
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for a 100 au distribution, plus one 'external' companion. 
Apart from the potential 'hot Jupiters', all these objects are 
on highly eccentric orbits. 

We found that, apart for the 'hot Jupiters', the proba- 
bility of ending up with a planet orbiting at a given distance 
from the central star increases with the distance. Thus, this 
scenario produces a decreasing number of planets as we go 
from 100 au down to 0.1 au. This is expected since a planet 
on a highly eccentric orbit has more chance of colliding with 
the central star as it gets closer to it. Paradoxically the ac- 
tion of the central star tends to clear out a large cavity 
around it. 

We now turn to the characteristics of some of the ex- 
trasolar planets detected so far. We restrict our attention 
to those with a projected mass Msini > 4.5 Mj, consis- 
tent with the larger masses expected for formation through 
fragmentation. There are 7 such isolated objects around 
HD 190228, UA3, HD 2222582, HD 10697, 70 Vir, HD 89744 
and HD 114762. All have an eccentricity larger than 0.12, 
and for 6 of them e > 0.33. The semi-major axis is in the 
range 0.88-2.5 for 5 of the objects, while the others have 
a — 0.3 and 0.43, respectively. The characteristics of these 
planets might be accounted for by our scenario though the 
two cases with small a would only occur with a small prob- 
ability. 

Amongst the 'hot Jupiters' detected so far, r Boo is 
particularly massive with M sin i ~ 4. This planet is clearly 
a candidate for being produced by a mechanism of the type 
we consider. 

The objects expelled as a result of the type of relax- 
ation process we consider may produce a population of freely 
floating planets ( Lucas & Roche 2000, Zapatero Osorio et 
al. 2000) which is several times larger than that of giant 
planets close to the central star. This population would be 
expected to be typically at least 10 times larger than the 
population of massive planets orbiting around the star, and 
depends on the initial number of planets in the distribu- 
tion. However, note that the population of planets orbiting 
central stars that went through a relaxation process of the 
type we consider may be significantly smaller than the to- 
tal if planets are also independently formed through a core 
accretion process. 

The model we have considered in this paper has of 
course some limitations. In particular, we have not included 
the effect of gas which is probably still significant in the early 
stages of evolution of these systems, and which may result in 
more planets orbiting at smaller distances from the central 
star. Also we have not considered the effects of Roche lobe 
overflow for the planets which have a close encounter with 
the central star (Trilling et al. 1998) which may significantly 
reduce their mass. 

The considerations above lead to the suggestion that 
there may be two populations amongst the extrasolar plan- 
ets detected so far. Some of the more massive objects may be 
produced through fragmentation of an envelope or a disc- 
like structure followed by dynamical relaxation. Other pre- 
dominantly lower mass objects could be produced in a disc 
as a result of the 'core accretion' model. We would of course 
expect to also have hybrid systems, in which both processes 
have occured. The relative importance of these processes 
may depend for instance on the physical parameters, such 
as the angular momentum content, of the parent cloud. 



If the scenario presented here is indeed effective in 
planet formation, we would expect additional massive plan- 
ets to be detected further from the central star, and an im- 
portant distribution of loosely bound, or 'free-floating' ob- 
jects associated with these systems. 
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RUN 1 




10 4 2xl0 4 3xl0 4 4xl0 4 5xl0 4 6xl0 4 




Figure 1. This figure shows the evolution of the semi-major axes (v/pper plot) and pericentre distances (lower plot) of the 5 planets in 
the system versus time (measured in units of Po) for run 1 in table hi The lines correspond to the different planets each having a mass 
5 X 10 — 3 M*. In this and other similar figures, a line terminates just prior to the escape of a planet. In this case R* = so that there 
was no tidal interaction with the central star. Note that the time resolution in this and similar figures is not fine enough for all close 
approaches to the star to be fully represented. 
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RUN 2 




10 4 2xl0 4 

time/P 

Figure 2. Same as Fig. U but for run 2 in table U. For this run N = 10 and R*/R max = 9.396 X 10 -5 . 
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RUN 3 
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Figure 4. This figure shows the evolution of the semi-major axis (solid line) and the pericenter distance (dotted line) of the innermost 
planet in run 3 vs. time (measured in units of Po) for a short interval after this planet enters a tidal interaction phase. This is a zoom 
on the curves displayed in figure ^. 
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RUN 5 
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RUN 6 
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Figure 7. The upper plot shows the evolution of the semi-major axes of the only 2 planets in the system which have not either been 
ejected or collided with the central star after about 180 Po versus time (measured in units of Pq) for run 6 in table |l| The lower plot shows 
the pericenter distance (dotted line) and the apocenter distance (solid line) for these planets. Here N = 100, R*/Rmax = 9.396 X 10 -5 
and the masses were selected uniformly at random in the interval (0; 10 — 2 M«). The innermost planet suffers a gravitational scattering 
by the more massive outermost planet when their orbits cross which results in an increase of its eccentricity and eventually a collision 
with the central star. 
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